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Abstract. We study the electrostatic Casimir effect and related phenomena in 
equilibrium statistical mechanics of classical (non-quantum) charged fluids. The 
prototype model consists of two identical dielectric slabs in empty space (the pure 
Casimir effect) or in the presence of an electrolyte between the slabs. In the latter case, 
it is generally believed that the long-ranged Casimir force due to thermal fluctuations 
in the slabs is screened by the electrolyte into some residual short-ranged force. The 
screening mechanism is based on a "separation hypothesis" : thermal fluctuations of 
the electrostatic field in the slabs can be treated separately from the pure image 
effects of the "inert" slabs on the electrolyte particles. In this paper, by using a 
phcnomenological approach under certain conditions, the separation hypothesis is 
shown to be valid. The phenomenology is tested on a microscopic model in which 
the conducting slabs and the electrolyte are modelled by symmetric Coulomb gases 
of point-like charges with different particle fugacities. The model is solved in the 
high-temperature Debye-Huckel limit (in two and three dimensions) and at the free 
fermion point of the Thirring representation of the two-dimensional Coulomb gas. The 
Debye-Huckel theory of a Coulomb gas between dielectric walls is also solved. 
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1. Introduction 

Two parallel metallic plates attract one another at zero temperature due to fluctuations 
of the quantum electromagnetic field in vacuum. This is the well-known Casimir effect; 
for a nice introduction see reference pQ , for an exhausting review see |2j . The extension 
of Casimir's result to arbitrary temperatures and to general dielectric plates was made 
by Lifshitz et al. with a subsequent treatment of delicate points by Schwinger et 
al. 4J. The Casimir force between two disconnected conductor objects of different 
shapes, like a rectilinear wall and a sphere, is also treatable |T| by using a method due 
to Derjaguin. The general form of the Casimir free energy for ideal-conductor walls 
of arbitrary smooth shapes was derived by Balian and Duplantier [3] using multiple 
scattering expansions (for a review, see [6 ). In the classical (or high-temperature) limit 
defined by the validity of the equipartitioning energy law, the Casimir forces between 
disconnected boundaries become purely entropic [7|, and as such depend only on the 
geometry of the boundaries. For accurate experimental measurements of the Casimir 
effect see jH]. 

There is a permanent interest in the Casimir effect and related phenomena in 
equilibrium statistical mechanics of classical charged fluids, i.e. models of non-quantum 
charged particles which do not incorporate the magnetic part of the Lorentz force due to 
charge currents. Although such purely electrostatic models ignore the magnetic degree 
of freedom (whose effect has the same magnitude as that of the Coulomb force in the 
case of conducting walls), they can be used in more complex experimental situations 
to reveal microscopic mechanisms behind macroscopic electrostatics. These classical 
electrostatic forces are also called van der Waals forces. The statistical models confined 
by walls can be divided, in general, into two sets: the semi-infinite systems, in which 
at least one of the spatial coordinates is unconstrained by the walls, and the fully- 
finite systems. The applied methods and observed phenomena usually depend on this 
classification. Here, we shall restrict ourselves to the semi-infinite confining geometries, 
namely the polarizable interface (figure 1) and the two-slab geometry (figure 2). The 
fully-finite geometries will be treated in a subsequent paper. 

The studied models include the possibility of the presence of an electrolyte between 
the conducting walls. It is generally believed that the long-ranged Casimir force between 
the fluctuating walls is screened by this electrolyte. Such an intuitively expected 
behaviour was supported longtime ago by an approximate Poisson-Boltzmann analysis of 
an electrolyte between two dielectric slabs jHj- A reasonable, but not rigorously justified, 
assumption was made: the electrostatic-field fluctuations in the walls, responsible for 
the long-ranged attractive Casimir force, can be treated separately from the pure image 
effects of the "inert" walls on the charged constituents of the electrolyte. Hereinafter, 
we refer to this assumption as the "separation hypothesis" . Considering then the pure 
effect of images on the electrolyte particles, a repulsive long-ranged force of the same 
amplitude as the Casimir one, but in the opposite direction, occurs in the electrolyte. 
Following the basic separation hypothesis, the two long-ranged forces cancel exactly 
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Figure 1. Polarizable interface. 



with one another; we shall refer to this cancellation phenomenon as the screening effect. 
By the use of the separation hypothesis, this screening effect was also observed for an 
electrolyte between two conducting walls . The presence of the repulsive long-ranged 
force in the electrolyte (modelled by the Coulomb gas) has been proven at arbitrary 
temperature, by using perfect-screening sum rules, for ideal-conductor jTT] and ideal- 
dielectric ^2] inert slab walls. Interestingly, in the special case when e = 1 (plain hard 
walls separating the electrolyte from vacuum), there is neither an attractive Casimir 
force between the slabs nor a repulsive force in the electrolyte induced by the particle 
images in the inert walls (in fact, there are no images). The corresponding field theory 
was developed in a series of works by Dean and Horgan j!3j . 

In references JTJJ [T^, and in a recent work ^B] about the microscopic origin of 
universality in Casimir forces, the conducting walls in empty space are modelled by 
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two slabs filled with a classical plasma of charged particles in thermal equilibrium. At 
large distances between the two slabs, the Casimir force (in the electrostatic regime) 
is retrieved as a result of perfect-screening sum rules. A natural continuation of the 
strategy is to replace the ideal conductor walls by a microscopic plasma system in more 
general physical situations. 

The aim of the present paper is to study Casimir and related effects in model 
systems with geometries presented in figures 1 and 2, in which the confining walls are 
made of a symmetric two-component plasma (Coulomb gas) of point-like particles with 
+/— unit charges in thermal equilibrium. The Coulomb gas is chosen for the sake of its 
simplicity; most of the obtained results can be easily generalized to mult i- component 
plasmas. In the presence of an electrolyte (which is simply another Coulomb gas with 
a different particle fugacity) between the walls, the microscopic description of the walls 
enables us to mimic coherently the effect of electrostatic fluctuations inside the walls 
and the image forces acting on the electrolyte particles, without any ad-hoc separation 
ansatz. The wall-plasma is not required to be ideal, there might be a finite nonzero 
correlation length which describes the correlation effects among the charged particles. 

If there is no electrolyte in the region between the two slabs (figure 2), the Casimir 
force is obtained in the limit of a large wall separation, together with a leading correction 
term due to the non-ideality of the conductor wall. If there is an electrolyte between 
the two slabs, the Casimir force is shown to be screened by this electrolyte, which 
confirms the validity of the separation (fluctuations-images) hypothesis. Interestingly, 
the residual short-ranged force between the slabs is always attractive. As a by- 
product of the treatment, we consider also the ideal-conductor regime of the walls when 
the correlation length of the charged particles forming the walls goes to zero. The 
surface tension of the electrolyte turns out to be independent of the electrostatic field 
fluctuations in the ideal-conductor walls. The analogous independence of the density 
profile and particle correlation functions was previously noted jTZl HE] • 

The spatial dimension of the studied models will be v = 2, 3. In dimension is, the 
Coulomb potential v at a spatial position r, induced by a unit charge at the origin 0, is 
the solution of the Poisson equation 



where s u = 2n u / 2 /T(is/2) is the surface area of the //-dimensional unit sphere. Explicitly, 



Here, r = |r| and a is a free length scale. The definition (jl.lj) of the Coulomb potential 
implies in the Fourier space the characteristic small- behaviour -O(k) oc 1/k 2 in any 
dimension. This maintains many generic properties (like screening JH1) of "real" 3D 
Coulomb systems. The interaction energy of charged particles {i,qi = ±1}, immersed 
in a homogeneous medium of dielectric constant = 1, is J2i<j QiQj^i^i — r j\)- 

The models are treated in the grand canonical ensemble characterized by the 
inverse temperature f3 and by the couple of equivalent (there is no external electrostatic 






otherwise. 
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potential) particle fugacities z+(r) = Z-(r) = z(r). The whole domain A on which 
the system is defined can be separated into disjunct subdomains, A = U a A^ a \ Within 
the grand canonical formalism, each subdomain is characterized by a constant fugacity, 
[r) = z^ a ' for r G A^ a \ which may vary from one region to the other. The choice 



z^ = corresponds to a vacuum subdomain with no particles allowed to occupy the 
space. The grand partition function is defined by 

oo i . N 



7V + ,iV_=0 +' 



;i.3) 



where N + (AL) is the number of positively (negatively) charged particles and N = 
N + + AL. The averaged particle densities are generated from H in a standard way as 
functional derivatives with respect to the fugacity field. At the one-particle level, one 
introduces the number density of particles of one sign, 

1 5E 



„„(r) (J25 q , qi 6(r-r i )) = --^-. (1.4) 



'5 V 

Due to the charge symmetry, n+(r) = n_(r) = n(r)/2 where n is the total density of 
particles. At the two-particle level, one introduces the two-body densities 

7v(r,r') = ($]<W( r - r^^r' - Tj) ) 

1 S 2 E 



;i.5) 



S5^(r)5zg/(r')' 

It is useful to consider also the (truncated) pair correlation functions 

n 3 (r)n g /(r') 

For the case of point-like particles, the singularity of the Coulomb potential ()1.2|) at 
the origin often prevents the thermodynamic stability against the collapse of positive- 
negative pairs of charges: in 2D for small enough temperatures, in 3D for any 
temperature. In particular, in 2D, the Boltzmann factor of a pair of +/— unit charges, 
r _/3 , is integrable at short distances provided that j3 < 2. To cross the collapse 
point (3 = 2, the pure Coulomb potential has to be regularized by a short-distance 
repulsion, e.g., a hard-core potential. In 3D, the short- distance regularization is needed 
at any temperature (except in the high-temperature limit described by the Debye- 
Hiickel approximation which gives a thermodynamically stable description for point-like 
charges) . 

Two methods are applied. We first use the Debye-Hiickel method of the Coulomb- 
bond chain resummation within the density Mayer expansion, valid in the high- 
temperature limit and applicable to both 2D and 3D systems. The standard approach 
dealing with inhomogeneous density profiles fK)\ |2*T] has to be supplemented by the first 
equation of the BGY hierarchy. We rather use the version developed in |2*2^ 12*3*] which 
does not require some additional information; we also introduce some technicalities 
which substantially simplify the calculation of the Casimir forces. We also use the 
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Debye-Hiickel method in the case of dielectric walls, using the separation hypothesis 
which does not require any microscopic model for these walls. The second method 
applies exclusively to the 2D Coulomb gas whose exact solution (Ursell functions) can 
be obtained just at the collapse temperature (3 = 2 by a mapping onto the Thirring 
model at the free-fermion point ^3 El]- The free-fermion model is solvable also in 
many inhomogeneous situations and the corresponding method can be adapted to the 
present problems. The fact that, in 2D, one can verify whether some features of the 
mean-field behaviour predicted by the Debye-Hiickel analysis persist also at a specific 
finite temperature is of importance. We add that, based on a recent progress in the 
integrable 1 + 1 dimensional Quantum Field Theory [23 EE], the thermodynamics of the 
2D Coulomb gas was solved exactly in the whole stability region (3 < 2, in the bulk [27] 
as well as for interface geometries with an ideal-conductor [2Hj and an ideal-dielectric 
[2*H] walls; however, for (3 < 2, the exact n-body densities are not known, and we shall 
not consider the case (3 < 2 in the present paper. 

The paper is organized as follows. In section 2, for the two-slab geometry (figure 2) 
with general dielectric walls, by using a phenomenological approach, it is shown, under 
certain conditions, that the electric field fluctuations can be decoupled from the image- 
particle forces. The separation hypothesis is then tested on the microscopic models, first 
in the Debye-Hiickel limit (section 3) and then at the free-fermion point of the Thirring 
representation of the 2D Coulomb gas (section 4). A brief recapitulation and concluding 
remarks are given in section 5. 

2. Phenomenology 

2.1. The separation hypothesis 

A general argument in favour of the separation hypothesis has been given in the case 
of conducting walls ^B]- After a minor modification, the same argument can be used 
for dielectric walls of arbitrary dielectric constant e. For the sake of completeness, we 
repeat the argument. 

We consider, in v dimensions, a Coulomb gas partially or totally bounded by walls 
W made of a dielectric with a dielectric constant e. We want to take into account the 
thermal fluctuations of the walls. It happens that a detailed knowledge of the internal 
structure of the walls is not necessary; we only assume that this internal structure is 
described by a set of coordinates {R}. Classical (non-quantum) statistical mechanics is 
supposed to be applicable. 

We shall need a correlation function defined for the wall system supposed to be alone 
in space. Let $w( r ) be the electric potential created by the walls at some point r and 
let (<lv( r ) < ^V( r '))o be a truncated statistical average computed with the Boltzmann 
weight of the walls alone. This correlation function, when both r and r' are outside the 
walls, can be computed by linear response theory, by a simple adaptation of reference 
[3*U] . Let us put a point charge q at r'. Its interaction with the walls results into the 
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addition to the Hamiltonian of a term 5H = q§w(r'). The average electric potential at 
r is changed by qG(r, r') with G determined by macroscopic electrostatics. For instance, 
in the case of one plane wall, G(r, r') = v(r — r') + [(1 — e)/(l + e)]v (r — r'*) where v is 
the Coulomb interaction (1.2) and r'* is the image of r'. 

Let — G* be that part of G due to the walls (not to q itself): 

-G*(r,T f )=G(r,T f )-v(T-T f ). (2.1) 

If q is infinitesimal, linear response theory states that — qG* = — P{& w (t)5H)q , therefore 

/?<<Mr)<Mr')>? = G*(r,r'). (2.2) 

Furthermore, we can assume that the fluctuations of $vy(r) are Gaussian, because, in 
a macroscopic description with a constant e, the response of the walls to the external 
charge q is linear even for a non-infinitesimal charge (if the response of a spring to an 
applied force is linear, it means that the spring is harmonic, and its spontaneous thermal 
fluctuations will be Gaussian). 

We now come back to the full system Coulomb gas plus walls, the partition function 
of which will be considered. Let Ho({R}) be the potential energy of the walls, {r} the 
set of the particle coordinates of of the Coulomb gas. The total potential energy of the 
system (Coulomb gas plus W) is 

H({v}, {R}) = ]T q^ w (r h {R}) + ]T q iqj v(n - Tj) + H ({R}), (2.3) 

i i<j 

where and are the charge and coordinates of the ith particle of the Coulomb gas. 
Let dr the phase space element (after integration on the momenta) of the Coulomb gas. 
The partition function Z of the total system can be written as 



Z = Z J dr ^exp 



-/3$><Mr*,{R}) 



exp 



-/3^g*g^(ri - 



Kj 



(2.4) 



where (• • -)o means the average over {R} with the weight exp[— (3H ], and Z is the 
partition function of the walls alone. Since the fluctuations of & w are Gaussian and 
(<Mr)> = 0, 



exp 



-/35>*w( r <>{ R }) 



exp 



exp 



1,3 



Thus, the partition function can be factorized (separated) as 

% — Z c gZ , 



(2.5) 



(2.6) 



where Z e s is the partition function of the Coulomb gas computed with the effective 
Hamiltonian 



H eft = -2E4 2G *( r " r i)+EM3 G ( r " r 3)' 

i i<j 



(2.7) 
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H e R is the standard Hamiltonian of a Coulomb gas in presence of inert dielectric walls. 
For instance, in the case of one plane dielectric wall, the first term of H e g is the 
interaction of each particle with its own image, while the second term is the interaction 
of each particle with the other ones and the images of the other ones. 

Let A({r}) be any microscopic quantity which does not depend on {R}. Using the 
same kind of steps as for the calculation of the partition function Z, one can show that 
this average can be calculated with if eff only: 



In particular, the densities and many-body correlations in the Coulomb gas can be 
calculated with if e fr and therefore are unaffected by the fluctuations in the walls. 

On the contrary, the total partition function ()2.6|) does involve Z , and therefore the 
fluctuations in the walls enter the total free energy; if a wall is displaced, the variation 
of the free energy, and therefore the force acting on this wall, will be affected by the 
fluctuations in the walls. 

In the above, simple manipulations allow to replace the partition function Z e $ of 
the Coulomb gas by the corresponding grand partition function E eS also computed with 
the Hamiltonian H c g. Therefore, the grand potential of the total system is 



where f2 e ff is the grand potential of the Coulomb gas computed with the Hamiltonian 
H e ff and F is the free energy of the wall system alone in space. 

Furthermore, we have considered the case when the Coulomb gas is made of charged 
particles in vacuum. A more realistic model for an electrolyte takes into account the 
solvent by picturing it as a continuous dielectric medium with a dielectric constant 65. 
Then, in the above, e means the ratio ew/^s where now ey/ is the dielectric constant of 
the walls. 

In the limit e — > 00, the results about conducting walls [T^ are retrieved. 
2.2. Casimir force between two dielectric walls 

In the present subsection, we give a simple argument for computing the Casimir force 
between two parallel dielectric walls separated by vacuum. We work in v dimensions. 
The walls, perpendicular to the x axis, modelled as media with a dielectric constant e, 
occupy the regions x > d/2 and x < —d/2. The slab —d/2 < x < d/2 is empty. Each 
point r is defined by its Cartesian coordinates (x, r ) where r 1 - is the — dimensional 
vector perpendicular to the x axis. We shall compute the average xx component of the 
Maxwell stress tensor [31] in the empty region. 

Let G(r, r') be the electric potential at r when a unit charge is placed at the source 
point r' in the empty region. Macroscopic electrostatics gives for G the equations 




off 



(2.8) 



^ — ^eff + Fq 



(2.9) 




if -d/2 < x < d/2 , 

if x < -d/2 or x > d/2 , 



(2.10) 
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with the conditions that G and e(dG/dx) [with e = 1 if —d/2 < x < d/2] be continuous 
at x — —d/2 and x = d/2, and that G —>■ when x — > ±oo. The standard technique 
is to use the Fourier transform G in the r space, defined as in ()3.15jl . A derivation 
similar to the ones in section 3 gives for G*, defined by 1)2.1)) . the Fourier transform G* 



G* = — j- cosh[Z(s — x') 



+ cosh[/(x + x 



1 + e 



1 - e 
1 + e 



e 2id _ l 



Jd 



-Id 



(2.11) 



The xx component T xx of the average Maxwell stress tensor at r is, from (|2.2j) with 
(<D w (r))=0, 



T xx (r) 



2s v 



d 2 
dxdx' 



G*(r,v') 



(2.12) 



When G* is expressed in terms of its Fourier transform (j2.11j) . 1)2.12)1 becomes 



f3T x 



{2kY- 1 Jo 



dll 

s u-l 



v-l 



1 + e 
1-e 



-i -l 



2ld 



(2n) v ' 1 (2dy VI + e 
where $ is the Lerch function 



, u, 1 



oo n—l 



n=l 



n" 



(2.13) 



(2.14) 



—T xx is independent of r. It is the force per unit area, measured along the x-axis, 
which acts on the right wall. This force is attractive. It is related to the free energy of 
the wall system F introduced in section 2.1 by 

1 dF 



— T 



(2.15) 



A dd 

where A is the of each wall. 

If v = 3, (2.13) is in agreement with Lifshitz [3]. In the limit of conducting walls 
e — > oo, it reduces to f3T xx = ((3)/ (Snd 3 ), where £ is the Riemann zeta function. This 
is indeed the correct electrostatic part of the Casimir force between conducting walls in 
the classical limit 01 • 

If v — 2, in the limit of conducting walls, f3T xx = 7r/(24<i 2 ), in agreement with a 
previous result 



3. Debye-Hiickel theory 



The weak-coupling (high-temperature) limit of a Coulomb system is described by the 
Debye-Hiickel theory. 
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3.1. General formalism 

We first consider the bulk regime when the Coulomb system is defined in the whole 
space, A = R u , and the fugacities {z q } of particles with charges q G Q as well 
as the corresponding particle densities {n q } are constant in space. A consistent 
way of deriving the Debye-Hiickel theory is to start with the diagrammatic Mayer 
expansion in density. Summing up the chain diagrams leads to the renormalized 
bonds, and the cluster integrals converge in the renormalized format (see, for instance, 
[23[27J). This is equivalent to writing the Ornstein-Zernicke (OZ) equations for the pair 
correlation function (jl.6|) with the direct correlation function replaced by —(3 times the 
corresponding interaction potential: 

h qiq2 (ri,r 2 ) = -Pq\<hv(\*i ~ r 2|) 

+ J2 / d u r 3 [-Pq 1 q 3 v(\r 1 - r 3 \)]n q3 (r 3 )h q3q2 (r 3 ,v 2 ) (3.1) 

The solution is of the form 

h qiq2 (r 1 ,r 2 ) = -Pqiq 2 G(ri,r 2 ), (3.2) 
where G obeys the integral equation 

G(r 1 ,r 2 ) = i;(|r 1 -r 2 |)-— / d"r 3 v(\ ri - r 3 |)G(r 3 , r 2 ). (3.3) 

S u J A 

Here, k is the inverse Debye length defined by k 2 = s u PJ2 q n qQ 2 - The integral equation 
(13. 3 j) can be transformed into a differential equation by taking the Laplacian with respect 
to ri . In this way one arrives at the usual Debye-Hiickel equation 

Al - k 2 ] G(n, r 2 ) = -8j(n - r 2 ) (3.4) 

which has to be supplemented by the vanishing boundary conditions for the screened 
Coulomb potential G when |r x — r 2 | — > oo. The bulk density-fugacity relationship can 
be obtained by using the formula 

f/V 1 
n q = z q exp I — lim \-G(r) + v(r)\ \ . (3.5) 

The specific grand partition function can be calculated from the standard relation 

d InS 

U « = Z «dz~ q |A| (3 ' 6) 
with the condition InH = at {z q = 0}. The bulk pressure p is given by ftp — In5/|A|. 
In the 2D Coulomb gas, the screened Coulomb potential is found in the form 

G{r) = K {Kr), k = {27T(3n) 1 / 2 [v = 2), (3.7) 
where K is a modified Bessel function. Using the small- a; expansion of K (x) |H2] 

K (x) = -In (|) - C + 0{x 2 lnx) (3.8) 
with C being the Euler's constant, the density-fugacity relationship (13 .5j) reads 

~z~aW = It J exp t ■ (3 - 9) 
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The bulk pressure is given by ftp = (1 — ft/4)n. This equation of state turns out to be 
valid for any ft < 2. 

In the 3D Coulomb gas, the screened Coulomb potential is found in the form 

G(r) = eM ~ Kr \ K={ATi(3n) 1 ' 2 (v = 3) (3.10) 
The density-fugacity relationship then reads 



n 



z + V2^{ftz) 3/2 \ . (3.11) 

The bulk pressure is given by ftp = n — k, 3 / (247r). 

The above Debye-Hiickel theory can be straightforwardly extended to inhomoge- 
neous situations when the system domain A is separated into disjunct subdomains, 
A = U a A^ a \ Each subdomain is characterized by constant fugacities, z q (r) = z^ for 
r G A*"), and the corresponding bulk densities n q a \ At the lowest level, the densities in 
the OZ equations (|3.1|) are taken constant (equal to the bulk ones) within each of the 
subdomains, n q (r) = n q a "> for r G A^ a \ The solution of the OZ equations is again of the 
form ()3.2|) . where the screened Coulomb potential G obeys the integral equation 

G(r 1; r 2 ) = vfln - r 2 |) - — I , d^drx - r 3 |)G(r 3 , r 2 ). (3.12) 

- s„ Jam 



Here, n a is the inverse Debye length for the subdomain defined by k 

s ^En n i Q '? 2 ' Equivalently, 



A: - k 2 ] G(r 1; r 2 ) = -sJin - r 2 ), n G A^, (3.13) 

where the spatial position of the source point r 2 is arbitrary. Since, according to (|3.12|) . 
G is proportional to v, these equations must be supplemented by the usual electrostatic 
conditions at each subdomain boundary dA^: G and its normal derivative with respect 
to the boundary surface d±G are continuous at dA^ a ' . The leading /^-correction to the 
particle densities in a given subdomain is determined by the inhomogeneous counterpart 
of equation ()3.5|) . 



n W(r)=4 a) 



exp j££ lim [-G(r, r') + v(r, r')]| (3.14) 



where r G A^. 



3.2. Polarizable interface 

The model is an infinite Coulomb gas of point unit charges q = ±1, separated by an 
interface (perpendicular to the x axis and localized at x = 0) into two half spaces x > 
and x < (see figure 1). Each point r G R u is thus defined by Cartesian coordinates 
r = (x, r ), where r 1 - is the {y — 1) dimensional vector component perpendicular to the 
x axis. The particle fugacities are, in general, different on each side of the interface, 
z q (r) = z when x > and z q {r) = zo when x < 0. The corresponding bulk (total 
number) densities and inverse Debye lengths will be denoted by (n, k) and (n , kq). Since 
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the system is translationally invariant in the perpendicular r^-subspace, the standard 
technique is to use the Fourier transform 

G(r 1; r 2 )= / ^L^^^explil-^-r^-)} (3.15) 

We want to solve equations ()3.13|) when, for instance, the source point r 2 = (x 2 , r^) 
has the coordinate x 2 > 0. In terms of the Fourier component G one obtains ordinary 
differential equations in one variable: 

/ d 2 



ydx 2 
dx\ 



k -/ \G(x 1 ,x 2 ;l) = —s„5(xi — x 2 ) 



/to - / 2 G(xi,x 2 ; 1) = 0, 



(xi > 0), 
Oi < 0). 



(3.16a) 
(3.166) 



We look for a solution, vanishing at X\ — > ±oo, of the form 
exp (— V K 2 + l 2 \xi — X 2 



G 



2Vk 2 + 1 2 



+Ai exp -Vk 2 + l 2 {xi + x 2 ) 



G = B[ exp 



/->•,- 



+ l 2 X\ — V K 2 + l 2 X 2 



(xi > 0), 
(xx < 0). 



(3.17a) 
(3.176) 



The constants Ai and Bi are determined by the conditions that G and dG/dxi are 
continuous at x\ = 0. The result is 



Ai 
B, 



Su ^pTp~-Jk 2 + p 



1 



(3.18a) 
(3.186) 

V^+P + sJk 2 + I 2 

The leading /3-correction to the constant particle density can now be calculated for 
x > by using a linearization in (|3.14j) . Considering x%,x 2 > 0, the first term on the rhs 
of (j3.17aj) contributes to the bulk density-fugacity relationship, while the second term 
implies 



?? 



n(x] 

More explicitly, one has 



X 



n\x) — n 



°° dll 



u-2 



kI + P- Vk^+I 2 



exp (— 2\/ k 2 + l : 



X 



(3.19) 



(3.20) 



for the density profile at x > 0. The same method applies when the source point 
i"2 = (x 2 , r^) has the component x 2 < 0, with the final result 



n (x) - n 



4(2< 



u-l 



K 2 + I 2 ^Tp + Jk 2 + P 



exp ( 2\l Kq + Px 



(3.21) 



for the density profile at x < 0. 
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For the present geometry, the grand potential Q - 
part and a surface part, 

^ = ~ ^ + + I^ItO, mo), 

where |<9A| stands for the interface area and 7 for the surface tension. The total number 
of particles is given by 



P In H is the sum of a volume 



(3.22) 



iv — z h z — . 

OZ OZq 

The surface parts of this relation read ^7] 



dz 



dx [n(x) — n] 



dz Q 



(3.23) 

(3.24a) 
(3.246) 



dx [n (x) — no] • 

The surface tension is zero in the homogeneous bulk regime, 

^j(z, zq) = when z = z . (3.25) 

At the lowest order in /3, zd z and z d zo can be replaced respectively by nd n and n d no . 
Having at one's disposal the density profiles, equations (|3.20p and (J3.21J) . after some 
algebra one finally arrives at 

h 



Su-l 



4(2tt; 



u-l 



dl l u - 2 ln 



+ p + 
aJk 2 + IWk 2 + P 



(3.26) 



Note that the argument of the logarithm behaves like 1 + 0(1 /I 4 ) for large /, so the 
integral is finite in both v = 2 and 3 dimensions of interest. When there are no particles 
at the half-space x < 0, i.e. k = 0, the relations 



(vac) 



K) 




32tt 



2 
3 



(3.27) 



define the surface tension of the plasma in contact with a vacuum (plain hard wall with 
e= 1) |33|. 

We now aim at studying the formula (J3.26|) in the ideal-conductor limit when the 
correlation length of the charged particles, say in the half-space x < 0, goes to zero, i.e. 
k,q — > 00. The analysis depends on the spatial dimension. 



In 2D, the large-K expansion of (j3.26j) results in a decomposition 



/37=^(4-vr)-^+/3 7ll 

07T O 



where 



a K 
/57mix = — , 

2tt Jo 



dk In 



VTTW + k 



(3.28) 



(3.29) 
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On the rhs of (|3.28|) . the first term corresponds to the vacuum surface tension of 
the Coulomb gas in the half-space x < (owing to k q » k, the relatively sparse 
Coulomb gas localized at x > is in fact seen as vacuum) and the second term 
corresponds to the surface tension of the Coulomb gas in the half-space x > in 
the presence of an inert ideal-conductor wall |2Hj- The last term, defined by (J3.29)) . 
represents the "mixing" of the two surface tensions. It can be further simplified, 
via successive integrations by parts, to the form 



Pin 




K J 



where E is a complete elliptic integral. In the limit k/k,q — > 0, one has 



Pin 



h' 



1 K ( K \ / K 

In — +0 — 

47TK \K Q J \K 



(3.30) 



(3.31) 



so that the mixing surface tension vanishes in the ideal-conductor limit. This means 
that the electrostatic field fluctuations in the wall do not modify the surface tension 
— k/8 calculated with an inert ideal-conductor wall. 

In 3D, the integral in equation ()3.26|) can be evaluated explicitly, 



16n 

•o, 



+ Kgln 



/ 2k 

\Kq + K 



+ K 2 \n 



2k 



K + K 



Pi 

In the limit k/kq 

/37=^(2.n2-l), 
where the mixing surface tension is given by the expansion formula 

Pi mix - 1 - 



K 



-In 



2k\ 

— )+Pln 
K 



12tt 





+ 




(-) 


~(-Y 


\k J 




\K J 



(3.32) 



(3.33) 



(3.34) 



Like in 2D, the first term on the rhs of ()3.33|) is the vacuum surface tension of the 
Coulomb gas in the half-space x < 0. The second term, the surface tension of the 
Coulomb gas in the half-space x > due to the conductor wall, diverges in the 
ideal-conductor limit. On the other hand, the mixing surface tension vanishes in 
the ideal-conductor limit. 



3.3. Two-slab geometry 

The model consists of two slabs (see figure 2) modelled by Coulomb gases with the same 
bulk characteristics (no, Kq). The slabs have thickness L and are separated by a distance 
al: 

A_ = {r = (x,r L )\x G X_ = (-d/2-L, -d/2),r L e R^ 1 } , 

A + = { r =(x,r 1 )|iGX + = (d/2, d/2 + L) , r x G R^ 1 } . 

There is an electrolyte, modelled simply by another Coulomb gas with different bulk 
characteristics (n, k), in the region (X = (-d/2, d/2)) x R"- 1 between the slabs; the 
choice n = k = corresponds to vacuum between the slabs. The remaining external 



Screening of classical Casimir forces by electrolytes 



15 



space to the slabs is supposed to be vacuum. The thickness L of the slabs is large enough 
to screen the boundary effects at x — id/2, so that the opposite boundaries have the 
density characteristics identical to those close to a plain hard wall. 

Let the source point r 2 = (x 2 ,r^) have the coordinate x 2 E X + . In terms of the 
Fourier component G, equations (J3.13)) take the form 

/ d 2 



\dx 2 
' cP_ 
dx\ 

dx\ 



- kI - l 2 j G{xi,x 2 ] 1) = -s v 8{xi - x 2 ) 

- K 2 -l 2 ^ G(xi,x 2 ;l) = 0, 

- - I 2 ) G(xi, x 2 ; 1) = 0, 



In the limit L —>■ oo, we look for a solution, vanishing at x\ 
G = 



2Jk 2 + 1 2 



exp -Jkq + Z 2 |xi - x 2 \ 



+Ai exp 



Kn + l 2 (xi + x 2 ) 



G = exp I — J Kq + l 2 x 2 I B[ exp ( v k 2 + l 2 x\ 



+Ci exp (— V k 2 + l 2 x\ 



G = Di exp 



Kq + l 2 (xi - x 2 ) 



{xt E X+), 

(x 1 E X ), 
(xi E XJ). 



(xx E X+), (3.35a) 

(xi E X Q ), (3.356) 

(xi E XJ). (3.35c) 
±oo, of the form 

(3.36a) 

(3.366) 
(3.36c) 



The constants Ai,Bi,Ci and D\ are determined by the conditions that G and dG/dxi 
are continuous at X\ = ±d/2. Introducing the auxiliary function 



W{1) = [(4 + I 2 ) + (k 2 + I 2 )] sinh (V^TPd 



+2 J 4 + IVk 2 + I 2 cosh ( Vk 2 + l 2 d 



one finds 



A 



K 



2W(/) J^TTi 



exp (\Jkq + l 2 d \ sinh (V k 2 + l 2 d 



B, 
C l 



2W{1) 
2W{1) 



l + l 2 + Vk^+P exp 



«g + 1 2 - v^T^ 



exp 



Kg + I 2 - VK 2 + I 2 



d 



(3.37) 

(3.38a) 

(3.386) 
(3.38c) 



Let us denote by f + the total average force per unit area acting on the slab 
A + as a whole. This force, directed along the x-axis perpendicular to the plates, is 
determined by the contact particle densities as follows. The electrolyte localized in 
the domain A contributes to this force by (3~ l n(d~ /2)e x , where e x is the unit vector 
along the x axis. The Coulomb gas localized in the considered slab A + contributes by 
p- 1 [-n Q (d+ j2) + n { ™ c) }e x , where is the particle density at x = d/2 + L, i.e. the 
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plain-hard- wall contact density for the Coulomb gas with the bulk density equal to n . 
The resulting force is given by 



PU 



'dr 



n 



n 



(vac) 

-no 



(3.39) 



The total force per unit surface acting on the slab A_, f_, is for symmetry reasons 
opposite to f + , f_ = — ff. The force f by which the slab A_ acts on the slab A + is 
determined as the difference of the force given by (j3.39|) and its asymptotic d — > oo 
separation value, 



Equivalently, 



where 



f(d) = f + (rf)-f + (oo) 



/3f = (Sn - 5n )e 2 



d 



dr 



8n = n\ — — n — 



5n c 



n TT 



d— too 



n ° IT. 



(3.40) 
(3.41) 

(3.42) 



d— >oo 



After some algebra one gets 
/3f (d) = 



2 u-l n (u-l)/2 T Jo 
1 



dl- 



X ■ 



W{1) 



exp 



Vk 2 + Pd) e x . 



(3.43) 



W(l) is defined by (J3.37|) . Note that the force between the slabs is always attractive. 
The large- distance asymptotic of this force depends on whether k = or k > 0. 

When there is vacuum between the slabs, i.e. k = 0, rescaling I into I = Kok, the 
amplitude / of the attractive force can be expressed as 



2 »-i n (»-m T ) Jo ( v ^pp + k 



dkk 



u-l 



exp (— kKod) 



Vl + k 2 + k) smh(knod) + 2kVl + k 2 exp(-kK d) 



-1 



(3.44) 



In the large- distance limit K d — > 00, one introduces the new integration variable 
k! G (0,oo) via k = k'/(n d), and then expands the integrated function in 1/(/c g(), 
with the result 



1 (v - i)r(i//2)c(i/) 



d" 



2 v ti v / 2 



1 



2z/ 

Kod 



O 



1 



(ko^) 2 



(3.45) 



The first term on the rhs of equation ()3.45|) corresponds to the electrostatic Casimir 
force in the ideal-conductor limit of the walls, the second term describes the leading 
(still long-ranged) non-ideality correction due to the finite correlation length 1/acq of the 
Coulomb gas forming the slabs. In the short-distance limit Kod — > 0, one finds from 
dUS) that 

M = , , f- Ak JZ,^ „ (3.46) 
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Explicitly, 

f Ko/(8tt), (i/ = 2), 

0/=^ , ; 1 ; (3.47) 

\«3/(247r), (i/ = 3). 

In the short- distance limit, the density of particles at x = d + /2 is in fact the bulk 
density n , so that (3f = n — Uq™^. According to [20], = n — Kq/(87t) in 2D and 



n, 



(vac) 



no — Kq/(247t) in 3D, which confirms the obtained results (|3.47j) . 



When there is an electrolyte between the slabs, i.e. k > 0, the amplitude / of 
the attractive force between the slabs, see equation ()3.43|) . can be expressed in the 
large-separation limit d — > oo as follows 

2(k 2 -k 2 ) 2 r°° dl I u ~ 2 Vk 2 + P ( r-a — - A 

(3f ~ ^ 7— TV / , -r exp ( -2V« 2 + l 2 d) . 3.48 

Introducing the new integration variable e (0, oo) via V^ 2 + ^ 2 — k + k/d and then 
taking the limit kg? — > oo, one finds the exponentially decaying force 

^/ - ^_ 1)/2 exp(-2«d). (3.49) 

If the slabs are made from ideal conductors (k — > oo) or they are empty (k = 0, i.e. 
plain hard walls with e = 1) this short-distance force takes the form 

Note that the parameter (2k) in the exponential decay corresponds to the inverse 
correlation length of the density-density distribution functions in the Debye-Hiickel 
limit [Slj. We conclude that the presence of an electrolyte between the conducting 
slabs induces screening of the long-ranged Casimir force. 

3.4- Plasma between dielectric walls 

We consider the case of a Coulomb gas between dielectric walls. We do not use any 
microscopic description of the walls and resort to the phenomenological separation (|2.9j) . 
The derivative dF / dd where F is the free energy of the wall system alone has already 
been calculated in section 2.2. We now use the Debye-Hiickel theory for computing the 
derivative dfl e s / dd where f2 e fj is the grand potential of the Coulomb gas computed with 
the Hamiltonian if e ff ; i>e. with inert dielectric walls. The geometry is the same as in 
section 2.2, except that now a Coulomb gas with the Debye inverse length k fills the 
slab -d/2 < x < d/2. 

We shall use (|3.14jl . thus we need the screened potential G. Let the source point 
r 2 = {^2, r 2") be in the Coulomb gas. The Fourier transform G of the screened potential 
is determined by the equations 
d 2 \ 

- k 2 -I 2 )G(xi,x 2 ]l) = sj(xi -x 2 ), (-d/2 < x x < d/2) , (3.51a) 




- I 2 } G(xi, x 2 ; 1) = 0, (x x < -d/2 or x x > d/2), (3.51&) 

oxf J 
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with the conditions that G and e(8G/dxi) (with e = 1 in the plasma region) be 
continuous at X\ = ±d/2 and that G —>■ when X\ — » ±oo. A derivation similar 
to the one in the previous section 3.3 gives, when both x\ and X2 are in the Coulomb 

gas, 



G(xi,x 2 ;l) = -^=J==|^exp (-Vk 2 + Z 2 |xi 



X 2 \ 



bi cosh 


\/k 2 + l 2 {xi + x 2 ) 


+ bf exp y—y/n 2 


+ l 2 dj cosh 


\/k 2 + / 2 (a;i - x 2 ) 


exp (Vk 2 + l 2 d 


- bf exp 




2 + l 2 d) 





where 



Vk 2 + l 2 -el 



(3.52) 



(3.53) 



v^TT 2 + el ' 

We now use (I3.52|) in ()3.14|) . The first term in the curly brackets of ()3.52|) 
corresponds to the bulk G; using the bulk relation ()3.5j) and linearizing the exponential 
in (j3.14|) gives 



/3n 



n{x) — n = — - 
Therefore 

rd/2 a r f \ i ? n 

ax[n[x) — n\ = 

d/2 2 



x 



d u - l l s u h cosh 


[2^K 2 + 1 2 


X 


+ 6 2 exp (-^K 2 + l 2 d) 


(27T)- 1 V« 2 + I 2 exp (vV + l 2 d) 


- 


-bfexp (-^K 2 + l 2 d) 





. (3.54) 



lu— 1 



(27T)"- 1 K 2 + I 2 



bi sinh KA 


2 + Z 2 dj 


+ bf^K 2 + l 2 dexp (-\/k 2 


+ Z 2 d) 


exp 


(Vrc 2 + l 2 d) 


- 6 2 exp (-v 7 ^ 2 + Z 2 d 





(3.55) 



The grand potential fi e g is the sum of a bulk part and a surface part (due to the 
walls) Qw This surface part is related to the surface part of the number of particles 
flXUni) by 

r d/2 

dx[n(x) — n] (3.56) 



d(-pn w ) 

n = A 



dn 



d/2 



since, at the lowest order in /3, we can replace zd/dz by nd/dn. 

In the limit d — > oo, Q w /A reduces to twice the surface tension 7, and ()3.56|) with 
the limiting form of ()3.55j) give 



n 



IV— 1 



1 h 



(3.57) 



dn 8 J (27t)^- 1 k 2 + / 2 ' 

The rescaling of I into I = nk makes apparent that the r.h.s. of ()3.57j) depends on n by 
a factor k u ~ 1 . Using nd/dn = (1/2)k<9/<9k, we can integrate (|3.57|) into 



Pi 



K 



IV— 1 



k 



1 VI + k 2 - ek 



4(1/ -1)7 (2ny- 1 1 + k 2 ^TTk 2 + ek' 



(3.58) 
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In 2D, performing the integral in ([3.58)1 gives 

2 , /l 



In 3D, for e = 1 



K 


71 


2^ 


~4 


K 


7T 


2tt 


~4 



VT^ 2 

2 



: tan 



: tanh 



-i 



1 + e 

e — 



1/2 



1\ 1/2 



e + 1 



if e < 1, 



if e > 1. 



(3.59) 



reproduces the z/ = 3 surface tension ([3.27)1 . For other values 
of e, the integral in ()3.58|) diverges. When e > 1, this divergence is to be expected for 
point particles, since the interaction between a particle and its image is attractive and 
produces a non-integrable Boltzmann factor. However, when e < 1, the divergence 
comes from the inappropriate linearization of the exponential in ([3.14)1 : a more careful 
treatment jHSl UBj gives a finite surface tension. 

For a finite distance d, the <i-dependent part of the grand potential per unit area 
is u) e ff(d) = {Vtw/A) — 27. After the substraction of the infinite-rf part of ([3.55)1 . using 
again nd/dn = (l/2)nd/dK, we now obtain 
d(M 1 



dn 



2(2tt) 



u-l 



X 



k[^(6 2 - \){k 2 + I 2 )- 1 + 2tfd{K 2 + I 2 )- 1 ' 2 } exp(-2d v / ^TF) 



(3.60) 



l-6 2 exp(-2cK/K 2 + / 2 , 
It can be checked that the numerator in the integrand of ()3.6()j) is the derivative with 
respect to k of the denominator. Thus, integrating ()3.60|) with the condition oj c q = 
when k = gives 
1 



/3w, 



cff 



2(2ny- 1 
-In 



1 



U In 



1 



1 



b] exp(-2dVK 2 + I 2 ) 



exj>(—2dl) 



(3.61) 



The second logarithm in equation (|3.6H) gives to the force —du e s/dd a long-range 
contribution opposite to the Casimir force in vacuum ([2.15)1 : the total long-range force 
associated to u = u cS + (F /A) vanishes. The cancellation (screening) effect does occur 
in the present model. Only a short-range (i-dependent force remains, associated to 

' ^ ]f-l7 1„ fl 1,2 1J./„2 1 ri\] (3 62) 



/3u 



2(2 



7T 



I In 



l-b 2 exp(-2dVn 2 + l 2 



Thus, in the 3D case v = 3, we have retrieved by our method a previously known result 

mm en]. 

In the limit nd — ► 00, the asymptotic behaviour of ([3.62)1 is 

1 



(3u 



J & v - l lb 2 exp{-2dVK 2 + l 2 ). 



(3.63) 



2(2tt)^- 1 

In this equation, only small values of l/n contribute to the integral. Thence, 
exp(— 2d\/ k 2 + I 2 ) ~ exp(— 2nd — dl 2 / k) and b 2 ~ 1. With these replacements in ()3.63j) . 

«("-i)/2 exp(_2«d) 



/3u 



2»(7Td)^ 1 )/ 2 



(3.64) 
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and the force per unit area has the asymptotic behaviour for large nd 
duo p- 1 k v exp(-2/td) 



(3.65) 



8d (47r/td)^- 1 )/ 2 ' 

This asymptotic behaviour has the remarkable feature that it does not depend on e. This 
was already apparent on (J3.49)) in the special cases e = 1 (kq = 0) and e — > oo (kq — > oo). 



4. The free-fermion point 

4-1- General formalism 

The 2D symmetric Coulomb gas is exactly solvable at the collapse coupling (3 = 2 by 
mapping onto the Thirring model at the free fermion point. We first shortly review the 
bulk regime. Both species of particles have the same rescaled fugacity m = 2naz [where 
a is the length scale introduced in equation 1)1.2)1 ] which has the dimension of an inverse 
length. The general formalism [TT) |2*4~] expresses the many-particle densities in terms 
of specific Green functions G qq '(r,r') (q,q' = ±). Because of the symmetry between 
positive and negative particles, one only needs G ++ and G_ + which are determined by 

(Ai - m 2 )G ++ (r u r 2 ) = -m5(n - r 2 ) (4.1) 

and 

G_ + (r 1; r 2 ) = -- ( A + G ++ ( ri , r 2 ). (4.2) 
m \oxi oyi J 

These equations have to be supplemented by the vanishing boundary conditions when 
| r i — r 2 | — > oo. In infinite space, the solution of equations 1)4.1)) and (|4.2|) reads 

m 

G ++ (ri, r 2 ) = —K (m\ri - r 2 |), (4.3a) 

G_+(ri,r 2 ) = - i i ATi(mri-r 2 , 4.35 

27r [ri — r 2 | 

where K and K\ are modified Bessel functions. The one-particle densities, given by 

n q = mG gg (r,r), (4.4) 

are infinite since Ko(mr) diverges logarithmically as r — > 0. This divergence can be 
suppressed by a short-distance cutoff R. We replace the point-like particles by small 
charged hard disks of diameter R and use a regularized form of ()4.4)1 for the total particle 
density n = n + + n_, 



2 ...2 , 2 



n = — K (mR) ~ — 

7T mfi^O 7T 



(4.5) 



Such a regularization is consistent in the sense that the perfect-screening sum rule is 
satisfied and the expected equation of state is reproduced [T7ll2^j . 

The above formalism can be generalized to inhomogeneous situations when the 
system domain A = U a A^ a \ Each subdomain A^ a ^ is characterized by a constant rescaled 
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fugacity, m(r) = m a for r G h^ a \ and the corresponding bulk density n a defined as a 
function of m a by (|4.5|) . The Green function G++ obeys the differential equation 

(Ai - m 2 a )G ++ (t 1 , r 2 ) = -m a 5(r a - r 2 ), n e A {q) , (4.6) 

where the spatial position of the source point r 2 is arbitrary. The Green function GL.+ 
is determined by C++ via equation (|4.2jl with m taken as the sub domain-dependent 
m(ri). The boundary conditions are that G ++ and G_+ must be continuous at each 
subdomain boundary dA^ a \ The one-particle densities are again given by (|4.4)1 . 

Polarizable interface 

In the geometry of figure 1, the rescaled particle fugacity is equal to m in the half-space 
a; > and to mo in the half-space x < 0. Let the source point r 2 = (x 2 , y 2 ) has the 
coordinate a; 2 > 0. In the Fourier representation ([3. 15)1 . equation ()4.6|) takes the form 

/ d 2 



\dx 2 
dx\ 



m —I \ G++(xi, x 2 ; I) = — m5(x\ — x 2 ), {x\ > 0), (4.7a) 



m 2 -/ 2 G++(si,x 2 ;0 = 0, (a* < 0). (4.76) 



We look for a solution, vanishing at x\ — > ±oo, of the form 

^++ = o i 7 , i exp(-v / m 2 + / 2 |xi -x 2 \) 
2ynr + / y 

+A 2 exp \-\lm 2 + / 2 (xi + a; 2 )l , (a* > 0), (4.8a) 



C++ = S, exp J ml + l 2 xi- Vm 2 + l 2 x 2 , (xi < 0) . (4.86) 



The constants A\ and -B; are determined by the conditions that G ++ and G_ + = 

[m(xi)] _1 (/ — d Xl )G ++ are continuous at x\ = 0. The final result is 

/ > ) 

\ m m / 

The particle density at x > is then given by 

fjl poo , . 

n ( x )= n+ — dl (A t + A A exp ( -2Vm 2 + l 2 x] , (4.10) 
n Jo v ' 

where n is the regularized bulk density related to the rescaled fugacity m via ()4.5j) . The 

same procedure applies when the source point r 2 = (x 2 ,?/ 2 ) has the component x 2 < 0. 

For symmetry reasons, one gets the particle density at x < in the form no (a;) = no+the 

last term of equation (|4.1()jl with the interchange m <-> mo- 

The surface tension 7 is again determined by the relations 

d(-(3j) f°° 

m — = / dx[n(x) — n\, (4.11a) 

dm 

dm 







= f° dx [n {x) - n ) , (4.116) 
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with the boundary condition /37(m,m ) = when m = m . The result is 

roo dl 



Pi 



In 



2^/m 2 + l 2 Jm 2 + I 2 



(4.12) 



lo 2tt \ymm + \/m 2 + l 2 ^jm\ + I 2 + I 2 ( 

When there are no particles in the half-space x < 0, i.e. mo = 0, ()4.12j) yields the 
surface tension of the Coulomb gas in contact with a vacuum, 

/3 7 (vac) ' 



TO 



m . 



(4.13) 



in full agreement with the result obtained in [T7] . 

In the ideal-conductor limit of the Coulomb gas in the half-space x < 0, the surface 
tension of equation ()4.12j) has the following large-m expansion 



mo 
Air 



2) + m 



2tt \m J v ; 



(4.14) 



The first term on the rhs of equation (|4.14j) is nothing but the vacuum surface tension 
[see (14.13)) ] of the Coulomb gas in the half-space x < 0. The second term, the surface 
tension of the Coulomb gas at x > induced by the ideal-conductor wall at x < 0, 
diverges due to the particle- image collapse at (3 — 2 [28] . 

4-3. Two-slab geometry 

The above method can be readily extended to the two-slab geometry of figure 2. We 
keep the notation of intervals X + = (d/2, oo), X = (—d/2, d/2) and X_ = (— oo, —d/2) 
from section 3.3. If the source point r 2 = (^2,2/2) has the coordinate X2 G X + , the 
solution of equations (|4.6J) is searched, in the Fourier space, in the form 

r - m ° 



2^/m 2 + 1 2 

+Ai exp 

exp 



exp ( — ymo + l 2 \xi — X2 

\Jml + l 2 (x 1 + x 2 ) 
ml + l 2 X2j B\ exp [\J m 2 + l 2 x\ 



(xi G X_ 



+Ci exp ^— V m 2 + l 2 X\ 



D t exp 



TOq + l 2 (xi - x 2 ) 



(Xl G Xq) 

(xi G X_ 



(4.15a) 

(4.156) 
(4.15c) 



The unknown constants A/, 1^, C/ and Z); are determined by the conditions that G ++ 
and [m(xi)}~ l {l — d Xl )G ++ are continuous at x\ = ±d/2. The particle density at x > d/2 
is determined by 

n (x) =n + ^J™dl (At + A- t ) exp (-2^ml + Px^j . (4.16) 

The same formalism applies when the source point r 2 = {x2,y 2 ) has the coordinate 
x 2 Glo, to obtain the density profile between the slabs. 

The force between the slabs is calculated by using the procedure outlined in section 
3.3, formulae ()3.39j) - ()3.43j) . Without going into technical details, we have found that 
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the force i(d), defined by (3.40), is always attractive. Its amplitude / as a function of 
the distance d between the slabs is given by 



md) = {m °~ m)2 fd/ ,w + ' 2 



* Jo mm + Vm 2 + l 2 y/m% + I 2 + I 2 

exp (-Vm 2 + l 2 d) , (4.17a) 



W(l) 
where 

W(l) = (mm + Vm 2 + l 2 ^fm 2 + I 2 + / 2 ) sinh (Vm 2 + l 2 d) 

+Vm 2 + l 2 \jm\ + I 2 exp (- Vm 2 + Z 2 d) . (4.176) 

The large-distance behaviour of f(d) depends on whether the rescaled fugacity m = 
or m > 0. 

When there is vacuum between the slabs, i.e. m = 0, the amplitude of the attractive 
force can be expressed in the limit m d — > oo as follows 

« = 55f{ 1 -pS55 + ((S5(p)}- (418) 

The first term on the rhs of equation 1)4.18)1 corresponds to the 2D electrostatic Casimir 
force in the ideal-conductor limit of the slabs, the second term corresponds to the 
correction due to the finite correlation length 1 / (2mo) of the Coulomb gas forming the 
slabs. It is interesting that this correction, evaluated at the (5 = 2 collapse point, has the 
same form (in terms of the correlation length) as the one evaluated in the Debye-Hiickel 
(3 — > limit, see formula (|3.45j) . 

When there is an electrolyte between the slabs, i.e. m > 0, the amplitude of the 
attractive force has the following large-ci behaviour 

(mo — m) 2 m 2 



d— »oo 



2 

m 



art ^3/2 ex P(- 2mrf ) form >0, 
exp(— 2md) for m = 



(4.19) 



. V Timd 

If the slabs are made from an ideal-conductor material with mo — * oo, this short-distance 
force is written as 

2 

m 

13 f ~ exp(-2mrf). (4.20) 

d->oo G\/'K{md) ii l' i 

These formulae confirm that, also at the special (3 = 2 coupling, an electrolyte between 
the ideal-conductor slabs causes screening of the long-range Casimir force. Comparing 
with the Debye-Hiickel result [see equations ()3.49j) and (|3.50|) ] one sees that the residual 
short-range force is temperature-dependent. Note that at (3 = 2 the residual force 
depends on the dielectric constant e of the slabs, since the results (14.19)1 for e = 1 
(mo = 0) and ()4.20j) for e — > oo are different. The non-dependence on e found in 1)3.50)1 
and ()3.65)) seems to hold only in the Debye-Hiickel regime. 
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5. Conclusion 

On fully microscopic models, in the classical (i.e. non quantum) regime, we have 
studied the properties of a Coulomb gas near a plane conducting wall, and two parallel 
conducting walls with or without a Coulomb gas between them. We have also considered 
the case of dielectric walls of arbitrary dielectric constant e, without a microscopic model 
for these walls, using macroscopic electrostatics and linear response theory. We have 
checked that the long-range Casimir force between parallel walls separated by a vacuum 
is screened into some (still attractive) residual short-range force when a Coulomb gas is 
present between the walls. 

All these calculations have been performed taking into account only the electrostatic 
part of the interactions. How to take into account the magnetic field for dealing, on a 
microscopic model in the classical limit, with the full particle-radiation interaction is an 
open problem. 
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